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Resonant inelastic X-ray scattering (RIXS) spectra of model copper- and nickel-based transition 
metal oxides are measured over a wide range of energies near the M-edge (hz/=60-80eV) to better 
understand the properties of resonant scattering involving shallow core levels. Standard multiplet 
RIXS calculations are found to deviate significantly from the observed spectra. However, by incor- 
porating the self-consistently calculated decay lifetime for each intermediate resonance state within 
a given resonance edge, we obtain dramatically improved agreement between data and theory. Our 
results suggest that these textured lifetime corrections are generic and significant for widely used 
resonant scattering spectroscopies. 



An accurate understanding of resonant X-ray interac- 
tions with matter is of central importance for current 
investigations into low energy many-body structure and 
dynamics with resonant elastic (REXS) and resonant in- 
elastic scattering (RIXS). Atomic multiplet calculations 
have provided the central basis for interpreting the so 
called "direct" scattering processes (non-shake-up pro- 
cesses) [iH^ by expanding upon a very precise formula- 
tion of intra-atomic electronic interactions [5]. Although 
these models have achieved remarkable success, current 
understanding of how they can be best implemented to 
study low energy material properties is limited by the 
complexity of many-body systems. 

In this Letter, we investigate RIXS spectra at the soft 
X-ray M2-M3 edges of SrCu02 and NiO, two relatively 
simple transition metal oxide systems that are well suited 
to analysis by atomic multiplet simulations. Our princi- 
ple finding is that the current continued use of just a 
single lifetime value to describe all of the intermediate 
resonance states at a given resonance edge is problematic 
in this energy range. We show that the atomic multiplet 
interactions can lead to more than 10 distinct lifetimes 
ranging from ~l-4eV for different intermediate states at 
the Ni M-edge. By incorporating self-consistent lifetime 
correction (SCLC) into the framework of existing mod- 
els, the variations in lifetime can be properly handled 
and the accuracy of resonant scattering simulations can 
be dramatically improved. Similar lifetime calculations 
have been performed to explain feature widths in mea- 
surements of core electronic states 0, 0] , however we will 
show that SCLC has more far reaching significance for 
the low energy degrees of freedom observed with RIXS 
and REXS. 

Nickel oxide is a model spin 1 antiferromagnet, with 
Ni at d8 valency. Low energy excitations in NiO M- 
edge RIXS spectra have been shown to correspond well to 
multiplet predictions SrCu02, with Cu at d9 and 

dlO valencies in the ground and intermediate states jioj . 
is an easily solvable model system from the standpoint 



of multiplet physics [11]. Comparison between these two 
materials not only provides an opportunity to test the 
validity of current atomic multiplet models in describing 
RIXS spectra involving extremely shallow core levels, but 
also offers a unique perspective into the detailed RIXS 
process at such low excitation energies. 

The need for intermediate state lifetimes to be self con- 
sistently calculated within numerical models can be seen 
in the stark contrast between M-edge NiO fluorescence 
yield spectra simulated with and without SCLC in the 
standard multiplet model (Fig. 1(b)). The difference be- 
tween corrected and uncorrected models becomes smaller 
in the harder X-ray regime, where lifetime is increasingly 
defined by core hole decay channels such as the L-edge 
3s-2p Coster-Kronig transition 12, [l3| that do not de- 
pend strongly on valence electron symmetries (Fig. 1(a)). 
Also, core hole shake-up excitations are relatively strong 
at the L-edge [1, [III, [H | , and have a nonlinear time depen- 
dence [15l4l7j that helps unify the lifetimes of multiplet 
RIXS features. Nonetheless, SCLC at the L2,3-edges can 
still improve the correspondence between theory and data 
for the measurement geometry in Fig. 1(c) by making the 
higher energy L2 peak broader and less prominent. 

M-edge measurements were performed at the beamline 
4.0.3 (MERLIN) RIXS endstation (MERIXS) at the Ad- 
vanced Light Source (ALS), Lawrence Berkeley National 
Laboratory. The data were recorded by a VLS based 
X-ray emission spectrograph equipped with a commer- 
cially available CCD detector [18]. Scattering intensity 
measured at each pixel was normalized to ensure that in- 
tensity in the final spectra accurately represents the den- 
sity of scattered photons per unit energy. Large cleaved 
single crystal samples were measured at a pressure of 
3xl0~^^ Torr at room temperature. The resolution- 
limited full width at half maximum of the elastic line 
is better than JE^25±2 meV for NiO M-edge measure- 
ments, and (5E'^200meV for the L-edge X-ray electron 
yield (XAS). Since the proximity of the [000] Bragg peak 
at the M-edge creates tremendous intensity at the elastic 
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FIG. 1: Shallow vs. deep core level resonant decay: (a) 

A diagram shows resonant excitation at the M- and L-edges 
of 3d transition metals, and the lowest energy dipole-allowed 
decay channel of intermediate states. Red lines indicate scat- 
tering processes that interact with electrons near the Fermi 
level, (b) Fluorescence yield measured at the M-edge of NiO 
is compared with (red) a standard multiplet calculation with 
inverse lifetime Fm^l.OSeV, and (blue) the same model using 
SCLC parameters oiCv{1.0eV) and Ccc=0.05eV. (c) (black) 
Total electron yield (representing X-ray absorption) measured 
at the L-edge is compared to (red) a standard multiplet cal- 
culation with constant intermediate state inverse lifetimes 
rL3=0.6eV and rL3=0.8eV, and (blue) the same calculation 
with SCLC valence transitions of Cv{0.3eV) and Ccc=0.2eV 
at L3. Ccc is increased to 0.5eV at L2 due to Coster-Kronig 
decay but Cv remains constant. 



line, making it challenging to measure fluorescence yield 
(FY) in a conventional way, FY in this study is mea- 
sured by integrating RIXS intensity within a 0.8-3.4eV 
energy loss window. Measurements at the M-edge of NiO 
were performed with in-plane polarization in the [001] 
scattering plane, and the photon beam had a grazing 
25^ angle of incidence to the cleaved [100] face. L-edge 
measurements were performed in the same configuration 
with a 45^ angle of photon incidence. All simulations 
are performed with photon geometries matching the ex- 
periment, and energy loss features are convoluted by a 
lOOmeV Lorentzian function. 

The direct scattering process of RIXS and REXS is 
described by the Kramers- Heisenberg equation as follows: 



Rf{E,hv,q) (X ^ 



E 



{f\Tl,\m){m\T^\g) 



hv — E„ 



(1) 



{E-EfY + {\Tff 
where g sums over the degenerate ground states, and m 



runs over intermediate states where a core hole is present. 
Incident photon energy is written as hv^ the excitation 
energy is and the inverse final state (|/ >) excitation 
lifetime is Fj. The photon perturbation Tq is given in 
the dipole approximation by a sum over spherical har- 
monics with = —1, 0, 1 for left^ linear {z) and right po- 
larized light. The sum over g and geometry-dependent 
sum over q' are required because the measurements are 
performed over a '^O.Smm long strip on the sample en- 
compassing many magnetic and structural domains fl9| . 
with no discrimination of the polarization of scattered 
photons. Unlike the L-edges, the M2 and M3 edges are 
merged together and not separately visible for late transi- 
tion metal oxides other than d9 systems such as cuprates 

mini [3. 

We have modified the equation from its most com- 
mon presentation by explicitly indexing the lifetime F^ 
of each intermediate state. Aside from this, our treat- 
ment follows that in Ref. [il,0-i,[il. Hartree-Fock val- 
ues are reduced to 80% and 70% for the inter-d-orbital 
and 3p-3d Slater-Condon parameters, and the d-orbital 
spin orbit coupling parameter is set to 0.083eV [20|. The 
intrinsic spin orbit splitting between 3pi/2 and 3p3/2 is 
set to 2.35eV, inflated by 30% from the nominal value 
of l.SeV for pure nickel due to reduced screening with 
positive 4s valence [ll|, [2l[ . All states are solved for by 
exact diagonalization on an atomic multiplet state basis. 
The e^ orbitals are modified with a crystal field pertur- 
bation of 10Dq=1.03eV and an external exchange field 
of J = O.leV to account for oxygen- mediated large spin 
antiferromagnetism. The exchange field is summed over 
all 12 [112] magnetic domain orientations. 

For the RIXS simulation in Fig. 2 (a, bottom), the in- 
verse lifetime F^ is set to a fixed value for all intermedi- 
ate states (F^=1.05eV), as in previous studies. Variants 
of this fixed-lifetime model have been used with consid- 
erable success at the leading edge of M-resonance [H, , 
and for the L-edge with distinct lifetime parameters F for 
L2 and L3 resonance [l|, [l, [l2|. When compared across 
the full incident energy range of M2/3 resonance however, 
the deviation from theory is striking. Although inelastic 
features labeled #1-4 are found at approximately correct 
energies, their predicted spectral intensities correspond 
very poorly with the experimental data in Fig. 2 (a, top). 

Integrated feature intensity is shown as a function of 
incident energy in Fig. 2(b-c), and is likely to be the best 
basis for rigorous comparison between the simulations 
and data, because the integration window is large enough 
that final states spread out on a ^O.leV scale by many- 
body effects will be fully counted. Features at higher in- 
cident energies are clearly too intense in the F^=1.05eV 
simulation and have lineshapes that are too sharp with 
respect to the data. Even at lower resonance energies 
near 66eV, the simulated intensity of features #2-3 is 
too strong relative to feature #1. To improve upon this 
modeling paradigm, we note that the inverse lifetime of 
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FIG. 2: Lifetime corrected simulation of 
M-edge scattering: (a) (top) Observed ex- 
citation peaks are numbered #1-4 on a plot of 
the incident energy dependence of M-edge RIXS 
scattering from NiO. A multiplet simulation is 
shown (middle) with SCLC and (bottom) with 
constant Fm^l-OSeV. (b-c) Spectral intensity 
is integrated from (dots) the RIXS data and 
(lines) simulations for excitation energy win- 
dows of (#1) 0.8-L3eV, (#2) L35-2.2eV, (#3) 
2.3-2.65eV, (#4) 2.75-3.25eV, (FY) 0.8-3.4eV. 
These correspond roughly to the 4 excitations 
predicted by multiplet theory, and to the in- 
elastic fluorescence yield. The experiment and 
simulations are scaled by the same factor in- 
dicated next to each curve. Statistical error 
from photon count and detector dark current 
is insignificant. 



an intermediate state is derived from the allowed decay 
channels, described to lowest order via Fermi's golden 
rule 



m] 



(2) 



where C is a constant. These matrix elements are ob- 
tained in the course of solving the Kramers-Heisenberg 
equation in multiplet RIXS models, to within a common 
multiplicative factor (Qj) that depends on the radial 
wavefunctions of electron shells between which a tran- 
sition occurs: 



(3) 



For NiO resonant scattering at the M-edge, fsd,sp in- 
dexes a final state in which the intermediate state has 
been terminated by the transition of a 3d electron into 
the 3p core level, and at the L-edge, /3s,2p would repre- 
sent a final state following from Koster-Cronig transition. 
Because decay channels that do not involve valence elec- 
trons give roughly constant contributions to the lifetime, 
it is convenient to reframe this as: 



r™ = Ccc+Yl (r-in - C'cc) I {fv 14 |m) I (4) 

q'Jv 

where Ccc is a constant representing the sum over 
transitions between core levels (and any other compet- 
ing decay channels), fy indexes all final states of va- 
lence electron decay, and the valence decay rate constant 
Cvi^min — Ccc) is defined such that the longest lived 
intermediate state has an inverse lifetime of F^^^. In the 



simulation of the L2 edge (Fig. 1(c)), using the values 
Cv{0.3eV) and Ccc=0-5eV causes the smallest interme- 
diate state lifetime to be F^i^=0.3+0.5=0.8eV. To incor- 
porate this state-by-state determination of intermediate 
state widths in our simulation for NiO M-edge RIXS, 
we use values of Cy(l.OeV) and Ccc=0.05. This as- 
sumes that competing decay channels not considered by 
the multiplet calculation, such as indirect RIXS excita- 
tions j23| and direct transitions from oxygen p-orbitals 
to the nickel 3p level, are very slow and constitute only 
O.OSeV out of the sum. 

Repeating the Kramers-Heisenberg multiplet RIXS 
calculation with intermediate state lifetimes determined 
in this self consistent way gives the SCLC spectrum 
shown in Fig. 2 (a, middle), which is in far better agree- 
ment with the experimental data. Curves tracing the 
intensity of inelastic features in the data and simulation 
show only minor discrepancies (Fig. 2(b)), and the calcu- 
lation very accurately reproduces the relative intensities 
of all features. The SCLC lifetimes do not just broaden 
the data, but also act as a physically necessary normal- 
ization factor of Equation (1). Without this normaliza- 
tion, the higher energy resonance around 70eV becomes 
far too intense, as well as being too sharp (Fig. 2(c)). 
Line shapes of resonance are generally broader at higher 
incident energies of hi^~70eV, and calculated interme- 
diate state lifetimes become shorter in that region. It is 
likely that the agreement between calculations and exper- 
iments can be further improved by including other ener- 
getic parameters such as the difference in nephelauxetic 
renormalization of Slater-Condon terms for e^ and t2g 
orbitals. 

Fig 3(a) illustrates how the diversity in M-edge inter- 
mediate state lifetimes comes about. For an isolated ion 
with no spin orbit interactions, few excitations are dipole 
allowed. However, intermediate state lifetimes are not 
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FIG. 3: The origin of M-edge resonant lifetimes: (a, 

top) Simulated SCLC RIXS scattering for an isolated Ni^^ 
ion with spin-orbit coupling disabled. Dipole excited inter- 
mediate state symmetries are highlighted, (a, middle) Simu- 
lated SCLC RIXS scattering from an isolated Ni^^ ion. (a, 
bottom) Simulated SCLC RIXS scattering from Ni^^ per- 
turbed by an octahedral crystal field. Insets show histograms 
of the intermediate state inverse lifetimes, (b-d) Histograms 
show intermediate state lifetime distributions calculated for 
cuprate, nickelate and cobaltate M2/3 RIXS scattering, us- 
ing the cross section parameters C3d,3p and Ccc fitted for 
NiO in all three cases. The inset of panel (b) shows that the 
M3 and M2 features of SrCu02 are accurately fitted with the 
predicted r=2.5eV lifetime. 



single valued (inset in Fig. 3 (a, top)). As an example, the 
^F intermediate state with angular momentum moment 
mL=+3 (i.e. L=3, mL=+3) is composed of a 3p core hole 
with mL=+l and a 3d hole with mL=+2. Selection rules 
dictate that the 3p core hole can be filled by d-electrons 
with mL=0 or +2, but there is a hole in the mL=-\-2 
3d-orbital, so only the mL=0 d-electron can transition 
to terminate the intermediate state. Thus, even with- 
out considering the detailed dipole matrix elements, one 
can see that the intermediate state will tend to take 
longer to decay than other intermediate states in which 
the p- and d-holes have opposed angular momentum. 

The addition of spin orbit coupling to the simulated 
ion causes new excitations to become allowed, but does 
not change the spectrum dramatically (Fig. 3 (a, mid- 
dle)). However, by breaking continuous rotational sym- 
metry with an octahedral crystal field (10Dq=1.03eV as 
for NiO), numerous non-zero transition matrix elements 



are created (Fig. 3 (a, bottom)). All features are shifted 
to higher energy by the crystal field perturbation, yet a 
remnant correspondence with the narrow and broad 
intermediate state lineshapes is visible in narrow fea- 
tures at ~66eV and broader features at ~70eV, respec- 
tively. This panel already closely resembles the experi- 
mental data, in spite of not including the antiferromag- 
netic exchange field. 

Another useful point of comparison is M-edge scatter- 
ing from SrCu02, which has full dlO intermediate state 
valence and thus only one intermediate state lifetime. Re- 
markably, using the parameters C3d,3p and Ccc derived 
from NiO results in an inverse lifetime of r^=2.48eV if 
starting valence is changed to d9 as for cuprates. This 
value correctly fits M-edge FY measured on SrCu02 (Fig. 
3(b, inset)), and is within error bars of the measured 
SrCu02 core hole lifetime value of r=2.5±0.3eV [11]]. 
The fact that the core hole lifetime of SrCu02 can also 
be accurately predicted from the multiplet calculation us- 
ing Equations (3-4) gives an extremely strong indication 
that core hole decay processes at the M-edge are accu- 
rately described by the model after inclusion of SCLC. 
Histograms of intermediate state lifetimes determined by 
multiplet calculation for NiO and CoO are presented in 
Fig. 3(c-d) showing that shifting the valence state away 
from full d-shell occupation increases the range of inter- 
mediate state lifetimes, and a lifetime-corrected model is 
likely to be even more important for 3d transition metal 
compounds with close to half-filled 3d shells. 

Resonant scattering involving shallow core levels is not 
as well studied as harder X-ray resonance, largely because 
of the weaker X-ray emission signal from strong Auger 
decay and poor penetration skin depths for soft X-rays. 
It is worth noting however that resonant self absorption 
is weaker relative to the non-resonant absorption back- 
ground at low photon energies, making close comparison 
between X-ray emission data and theory more practical 



111 . l24j . The corrections we identify are also relevant to 



harder X-ray regimes, and their effects in RIXS, REXS 
and FY scattering are enhanced by the need for accurate 
lifetime values to normalize the scattering equation (Eq. 
(1)) and identify the onset of core hole shake-up excita- 
tions 0, [l3 • Significant quantum interference between 
atomic multiplet scattering channels is also moderated 
by the intermediate state lifetimes [llj . These and other 
factors can generate tremendous complexity in theories 
of resonant scattering, thus it is essential to character- 
ize model systems such as M-edge NiO to identify the 
underlying scattering physics. 

In conclusion, through comparison of M-edge resonant 
inelastic X-ray scattering data on model systems NiO 
and SrCu02 with an atomic multiplet model, we have 
shown that intermediate state lifetimes for M-edge reso- 
nant inelastic scattering can be extremely diverse over a 
small incident energy range. By applying self consistent 
lifetime corrections (SCLC), the full inelastic resonance 
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profile of NiO is simulated with greatly enhanced accu- 
racy. Evaluating the case of CoO shows that the same 
principle can be expected to apply for resonant scattering 
on other compounds with shallow core levels, and deter- 
mination of these lifetime factors is critical for model- 
ing resonant soft X-ray probes of low energy many-body 
states. 
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